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Supernovae are rare nearby, but they are not rare in the Universe, and all past core-collapse super- 
novae contributed to the DifTuse Supernova Neutrino Background (DSNB), for which the near-term 
detection prospects are very good. The Super-Kamiokande limit on the DSNB electron antineutrino 
flux, cj){Ei, > 19.3 MeV) < 1.2 cm~^ s~^, is just above the range of recent theoretical predictions 
based on the measured star formation rate history. We show that the Sudbury Neutrino Observatory 
should be able to test the corresponding DSNB electron neutrino flux with a sensitivity as low as 
(f>(22.5 < Eu < 32.5 MeV) ~ 6 cm~^ s~^, improving the existing Mont Blanc limit by about three 
orders of magnitude. While conventional supernova models predict comparable electron neutrino 
and antineutrino fluxes, it is often considered that the first (and forward-directed) SN 1987A event 
in the Kamiokande-II detector should be attributed to electron-neutrino scattering with an elec- 
tron, which would require a substantially enhanced electron neutrino flux. We show that with the 
required enhancements in either the burst or thermal phase fluxes, the DSNB electron neutrino 
flux would generally be detectable in the Sudbury Neutrino Observatory. A direct experimental 
test could then resolve one of the enduring mysteries of SN 1987A: whether the first Kamiokande-11 
event reveals a serious misunderstanding of supernova physics, or was simply an unlikely statistical 
fluctuation. Thus the electron neutrino sensitivity of the Sudbury Neutrino Observatory is an im- 
portant complement to the electron antineutrino sensitivity of Super-Kamiokande in the quest to 
understand the DSNB. 

PACS numbers: 97.60.Bw, 98.70.Vc, 95.85.Ry, 14.60.Pq 



I. INTRODUCTION 

The core-cohapse death of a massive star and the sub- 
sequent type II optical supernova occurs at a rate of ~ 1 
per second in the Universe, each releasing a prodigious 
blast of ^ 10^® neutrinos and antineutrinos. While a 
burst of many neutrino events would be detected from a 
Milky Way supernova l), the expected occurrence rate 
is only ~ 3 per century 2] . For supernova as far away as 
10 Mpc, it should be possible to reliably detect neutrinos 
just one or two at a time, perhaps as often as once per 
year, with present and proposed detectors Q. However, 
most supernovae are vastly farther, such that the ex- 
pected number of neutrino events detected per supernova 
is <C 1. When weighted with the total supernova rate, 
however, the prospects are encouraging for the detection 
of the Diffuse Supernova Neutrino Background (DSNB). 
The Super-Kamiokande (SK) limit 4] on the flux of 
DSNB electron antineutrinos, (j){E^ > 19.3 MeV) < 1.2 
cm^^ s~^, isjust above the range of recent theoretical 
predictions 0, iZj 8]. In particular, the calculation 
based on the most recent astronomical datapredicts that 
the DSNB is on the verge of detectability [g • 

Searching for the DSNB signal requires the detection 
of an excess event rate over backgrounds, exploiting the 
different energy spectra of the signal and backgrounds. 
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Due to the near-isotropy of the scattered positrons in 
the reaction De + p ^ e'^ + n and the large distances to 
the supernovae, the detections of DSNB events will not 
be associated with particular optical bursts. If Super- 
Kamiokande is enhanced by the addition of gadolinium, 
as proposed by Beacom and Vagins, the detector back- 
grounds would be greatly reduced, allowing the clean de- 
tection of as many as 6 events per year j^, E| ■ Detec- 
tion of the DSNB may be the first detection of neutrinos 
from beyond 1 Mpc, and the second detection of super- 
nova neutrinos. In just a few years, the yield from SN 
1987A could be surpassed. This will provide valuable 
insight into supernova physics, neutrino properties, and 
the history of cosmological massive star formation. 

In the following, we discuss the prospects for detection 
of flavors of the DSNB other than f/g, which is the domi- 
nant yield in SK. Specifically, we propose using the Sud- 
bury Neutrino Observatory (SNO) to study the DSNB 
electron neutrino flux. SNO has a unique sensitivity to 
the electron neutrino flux though the charged-current in- 
teraction with deuterons, i^e + d —^ e~ +p+p. By combin- 
ing this charged-current reaction with neutrino-electron 
elastic scattering and the neutral-current breakup reac- 
tion 1/ + d I' + n + p, SNO has studied in detail the 
^B solar neutrino spectrum TT], with analysis of the hep 
spectrum forthcoming. Using the theoretical hep spec- 
trum and the measured atmospheric backgrounds above 
20 MeV at SK 111, we estimate the DSNB backgrounds 
at SNO. We then show that SNO could very likely reach 
a sensitivity as low as 0(22.5 < < 32.5 MeV) ~ 
6cm~^s~^. This would improve the existing limit from 
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Mont Blanc by about three orders of magnitude [l^l , and 
provide an important complement to the SK D^. limit. In 
addition to the importance of studying the emission 
for supernova models, it is also crucial for testing neu- 
trino mixing, as stressed by Lunardini and Smirnov |l3j| . 

There is an additional motivation to test the DSNB 
flux separately from v^^ stemming from the observed 
angular distribution of the SN 1987A data. Two well- 
studied features of this data that still stand out are that 
the first event in the Kamiokande-II detector was scat- 
tered forward, and the time-integrated angular distribu- 
tions in both the Kamiokande-II and 1MB detectors were 
both generally more forward and less isotropic than ex- 
pected from + P — > e+ -t- n events alone. Both features 
have typically been explained as either an increase in the 
neutrino-electron scattering rate beyond theoretical pre- 
dictions, or a statistical fluctuation ,14j. With the SN 
1987A data alone, this debate cannot be resolved. 

Increasing the Ue flux to the extent that these features 
of the SN 1987A data were a probable outcome is in con- 
flict with standard models of supernova neutrino emis- 
sion. However, given the fact that the models generally 
fail to produce successful explosions, the possibility of a 
serious misunderstanding of supernovae, perhaps due to 
new physics, cannot be excluded. Here we point out that 
such an increase is testable, in that it would generally 
lead to the DSNB flux being large enough for SNO to 
detect. The absence of a signal would greatly strengthen 
the case that these features of the angular distribution of 
the SN 1987A data were due to a statistical fluctuation. 
Resolution of this point is also important for the future 
interpretation of the SK DSNB results. 

In Section ^ we review the DSNB predictions and 
limits; in Section llTTI we calculate the capabilities of SNO 
to detect the DSNB Uf. flux; in Section llVl we apply these 
results to discuss the future constraints on the electron 
neutrino emission from SN1987A, and then we finish with 
our conclusions. 



II. PRESENT PERSPECTIVE ON THE DSNB 

A. Emission per Supernova 

Stars greater than ~ 8Mo are able to burn elements by 
successive nuclear fusion reactions until iron is reached, 
beyond which point no further energy generation is pos- 
sible Once the core has formed about a Chan- 
drasekhar mass of iron, not even electron degeneracy 
pressure can support it under the weight of the stellar 
envelope, and it collapses until it reaches the density of 
nuclear matter, bounces, and forms an outgoing shock. 
If energetic enough, the shock will eject the envelope and 
cause the optical supernova, leaving behind a neutron 
star; if not, the core and envelope will collapse into a 
black hole. In either case, the 3 x 10^^ erg of gravitational 
binding energy release is dominantly radiated away by all 
flavors of neutrinos and antineutrinos over about 10 s (if 



black hole formation occurs sooner than this, then the 
neutrino emission is less). It is usually assumed that this 
total energy release is shared more or less equally 
among the six flavors of neutrinos and antineutrinos. For 
recent updates and reviews on supernova neutrino emis- 
sion, see Refs. 0,^^. The spectroscopically classifled 
type II, lb, and Ic supernovae (hereafter SNII) result 
from massive star core collapse, and lead to the emission 
of ^ 10^^ neutrinos and antineutrinos. Type la super- 
novae (hereafter SNIa) result from the sudden ignition of 
a white dwarf accreting material from a companion, and 
are not accompanied by comparable neutrino emission. 

In core-collapse supernovae, neutrinos must diffuse out 
of the proto-neutron star, so that they should decouple 
with nearly thermal spectra, characteristic of their sur- 
face of last scattering. Since and and their antipar- 
ticles only have neutral-current interactions at these low 
energies, they should decouple at the smallest radius and 
largest temperature. Noting the charged-current interac- 
tions of Vf, and f'e, and that the medium is neutron-rich, 
one conventionally expected a hierarchy of temperatures 
like T^^ ~ 4 MeV, Tp^ ~ 5 MeV, and T^^ ~ 8 MeV. 
More recent work indicates both lower expected temper- 
atures, as well as lower differences between flavors |16j . 
Assuming equipartition of the total energy among the six 
flavors, the spectrum per flavor is approximately 
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Neutrino oscillations may mix the spectra of different fla- 
vors, and importantly, can increase the average energy of 
the more easily detectable and (for example, see 
Ref. ^3)' return to this point below. 



B. Integration over Past Supernovae 

It has been known for some time that the star forma- 
tion rate was larger in the past, and in particular, was 
about an order of magnitude larger at redshift 2 ~ 1 
than it is today, likely increasing more slowly at larger 
redshifts, and first turning on at an uncertain redshift 
2^5 (these data have been comprehensively reviewed 
by Hopkins In the past year, these results have 

been markedly improved by new data, especially from 
the GALEX ultraviolet satellite These and other 

recent data IE Els, while confirming the basic pic- 
ture above, also indicate that the absolute rates are on 
the high side of past estimates, due to larger (and better 
understood) corrections for obscuration by dust. 

Interestingly, this recent trend in astronomical analy- 
ses to favor larger dust corrections and larger star forma- 
tion rates can be very usefully constrained by the upper 
limits on the DSNB. In Ref. [g, Strigari ei al. developed 
a "Concordance Model" for the star formation and su- 
pernova rates which accounts for these and other data. 
This model is characterized by the star formation rate 
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per comoving cubic Mpc today, Rsfr{z = 0) 



0.02 



Mpc~^ yr~^, and the growth rate Rsfr{z) / Rsfr{z — 
0) = (1 + z)'^ , with f3 = 2, assumed to hold to 2; = 1, 
after which a constant rate was assumed. These choices 
of Rsfr{z — 0) and /3 are shghtly conservative, rela- 
tive to the GALEX data (see Fig. 1 of Ref. 8]). For 
larger z, the neutrinos are typically redshifted below 
the detector threshold, and so have relatively little ef- 
fect on the result. Taking into account the stellar ini- 
tial mass function, this converts to a core-collapse rate 
of Rsmiiz 0) = 2.5 x 10"'' Mpc'^ yr"!, with the 
same functional form, since SNII stellar lifetimes are 
short An important test of the Concordance Model 
is that it agrees with measurements of the SNII rate (and 
the SNIa rate too), confirming the choices of conversion 
factors, and requiring that the fraction of failed super- 
novae directly forming black holes must be relatively 
small 0. The neutrino emission parameters assumed 
were 3 x 10^"^ erg, equipartitioned among flavors, and 
normal-hierarchy mixing of T = 5 MeV z/g with T — 8 
MeV z/^ and Vt- For smaller temperatures and/or differ- 
ences between temperatures, the predicted flux would be 
somewhat less. 

Using the Concordance Model Q results for the SNII 
rate convolved with the neutrino emission per supernova 
dN/dE, the expected DSNB differential flux d(j)/dE is 

d<t>{E) f„ , ^dN{E{\ + z)) ,^ .dt^ 

where the integral runs from redshift zero to at least z ~ 
1, beyond which there is little contribution. To determine 
dt/dz, we use a A-CDM cosmology, = 1 ^ ^^A = 0.3, 
and Hq ^ 70 km s"^ Mpc"^, so that 
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Our results are in good general agreement with other 
recent calculations up to different choices of the in- 
puts d, IE ; these results and others are nicely reviewed 
by Ando and Sato The key input is the star forma- 
tion rate at redshift zero, and the most recent astronom- 
ical measurements used in Ref. Q are both larger and 
more precise than those used in Refs. H, IE S IM| ■ In 
Fig. n we show predictions for the DSNB i/g fluxes, for 
a few possible choices of the i/g temperature, without in- 
cluding any specific neutrino oscillation scenario yet. 



C. Existing Limits on the DSNB 



The SK upper limit on the DSNB flux of is by far 
the most stringent for any flavor, and is 1.2 cm~^ s~^ 
for energies above 19.3 MeV 0|. (With increasing en- 
ergy, the signal falls and the background rises, and so 
the flux at energies beyond about 30 MeV is less im- 
portant.) The detection channel in SK is inverse beta 
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decay on free proton targets, 
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and in 
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FIG. 1: Predictions for the DSNB spectra assuming 
T = 4, 6, or 8 MeV (solid curves), along with the approximate 
published Mont Blanc limit and the projected SNO sensitiv- 
ity (solid lines). The solar neutrino Ve spectra, which are 
irreducible backgrounds for this study, are also shown. The 
atmospheric neutrino fluxes are not shown here, but their ef- 
fects on the detectable spectra are shown in Fig. |5| 



the existing analysis, only the positron is detected 
At the present level of uncertainty and low statistics, the 
recoil-order corrections to the cross section and kinemat- 
ics [2^ can be ignored, so that E^ — Ei^ — 1.3 MeV, 
where Ee is the positron total energy, following the SK 
convention. Other singles rates in the same energy range 
as the signal create serious backgrounds. Below Ee — 18 
MeV, uncut muon-induced beta radioactivities and solar 
neutrino events overwhelm the DSNB signal. At higher 
energies, the dominant background is from the decays of 
sub-Cerenkov muons, produced by atmospheric and 
j>p interacting inside the fiducial volume. These non- 
relativistic muons quickly lose energy and decay at rest, 
with the relativistic electrons and positrons following the 
well-known Michel spectrum. Since their spectrum is 
known, and the total rate can be normalized from the 
data, the results are nearly independent of atmospheric 
neutrino flux uncertainties. The atmospheric Ve and De 
fluxes create a smaller background that can be well fit at 
energies above the Michel peak. Figure 2 of Ref. 0| is 
a clear illustration of these backgrounds and their large 
impact on the DSNB signal sensitivity. Importantly, re- 
cent DSNB calculations |E IE 011, 13 are not far below 
the SK limit, and the Concordance Model result based 
on the latest astronomical data is especially close 8] . 

The prospects for the future may be even brighter. 
Beacom and Vagins proposed that if SK were en- 
hanced by the addition of dissolved gadolinium trichlo- 
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ride (GdCla), then it would be possible to tag the neu- 
trons produced in the signal reaction, greatly reducing 
the background rates @- Gadolinium has a huge cross 
section for radiative neutron capture, and the produced 
gamma rays would Compton scatter electrons, produc- 
ing detectable Cerenkov light. With much-reduced back- 
grounds, a more favorable energy range could be used 
(beginning near 10 MeV), and SK could cleanly detect as 
many as 6 DSNB events per year; see Fig. 1 of Ref. 0. 
This could lead to the first detection of DSNB neutrinos. 
While gadolinium compounds have been commonly used 
in oil-based detectors, Ref. 9] was the first to propose 
that dissolved gadolinium could be added to a water- 
based detector, and that this is the only cost-feasible op- 
tion to develop a neutron-sensitive detector as large as 
SK, or ultimately at the 1-Mton scale ^2,]. The research 
and development work continues to be very encourag- 
ing, supporting the claims that the required standards 
of solubility, ease of use, water transparency, radiopu- 
rity, effects on detector materials, cost, and safety will 
be met 0,0. As a full system test, in late 2005 Vagins 
will add GdCla to the 1-kton water-Cerenkov detector at 
KEK (a scale model of SK, and a former near detector 
of the K2K long-baseline neutrino experiment [2^). 

The only published limits on the DSNB fluxes of other 
flavors come from the Mont Blanc experiment 0| . Their 
limit on v^^ flux is 0(25 < < 50 MeV) < 6.8 x 10^ 
cm^^ s^^, and is based on i/g charged-current interactions 
with ^^C. Since the DSNB spectra are quickly falling, 
even when weighted with the detection cross section, 
nearly all of their sensitivity would have come from the 
beginning of the above energy range. To make a more di- 
rect comparison with the sensitivity we derive for SNO, 
we simply assume that the Mont Blanc limit also applies 
to an interval of width 10 MeV, i.e., from 25 to 35 MeV. 
In Fig. n we show both the Mont Blanc limit and the 
projected SNO sensitivity as if they were constant over 
10 MeV, which is an approximation (e.g., note Fig.O. 

There is also a Mont Blanc limit on the i/^ and v-^ 
flux, based on neutral-current interactions with ^^C, and 
this is 0(20 < < 100 MeV) < 3 x 10^ cm'^ s-\ 
with a similar limit for v^^ and Vt. These are very weak 
limits (for comparison, the solar ^B flux is 5 x 10^ cm~^ 
s~^), and could be greatly improved. For example, a 
DSNB flux as large as the Mont Blanc limit would cause 
a huge excess of neutral-current deuteron breakup events 
in SNO, at least 10^ neutrons per day. It may be sufficient 
to set limits on the fluxes of DSNB and ly^ and use 
knowledge of the neutrino mixing angles to deduce limits 
on the DSNB fluxes of i'^, I'r, and their antiparticles. 

III. DSNB DETECTION IN SNO 

A. Detection of + d + p + p 

SNO, the flrst water-Cerenkov detector to use 
deuterons as a target for astrophysical neutrinos, is based 



on 1 kton of heavy water, D2O [27I l2^. The deuterons 
are targets for the charged-current reactions Vf, + d 
e~+p+p and De + d e++ri + n; electrons and positrons 
are detected by their Cerenkov light, as in light-water de- 
tectors like SK. Additionally, SNO has the ability to de- 
tect neutrons, so that it can measure the neutral-current 
reactions v + d v + p + n and i? -\- d D -\- p + n, sen- 
sitive to all flavors of neutrinos and antineutrinos, and 
separate them from the two charged-current reactions. 
Three neutron detection techniques have been used in 
SNO: radiative neutron capture on deuterons in the pure 
D2O phase, radiative neutron capture on chlorine in the 
D2O plus dissolved NaCl phase, and capture on discrete 
•^He-based neutron counters at present [S^, H^- In all 
cases, neutrons are simply counted, with no energy or 
direction information available. 

In a search for v^, SNO has the advantage of the large 
and spectral (good fidelity between incoming neutrino 
and electron energy) cross section on deuterons; SK de- 
pends on the much less favorable neutrino-electron scat- 
tering cross section. In contrast, in a search for Dg, SK 
has the advantage of a much larger detector size, as well 
as a more favorable cross section on free protons. The 
other key features of SNO are the existence of a good 
neutral-current detection channel and the ability to de- 
tect neutrons. While the former is unique to D2O, SK 
may soon have the latter ability |^ 0| , leveraged by its 
much greater fiducial mass of 22.5 kton. 

Other authors have considered the SNO sensitivity to 
the DSNB D^, noting that the very clean signal coinci- 
dence of a positron and two separate neutrons should 
allow better background rejection than in SK pol 
With limited exposure and low neutron detection effi- 
ciency (in the pure D2O phase), the SNO results so far 
are not very restrictive, but significant improvements are 
expected |31|. 

What we are proposing here, which has not been noted 
before, is that SNO should be able to make a unique con- 
tribution by exploiting its sensitivity to DSNB i/g. Since 
this will be just a singles search (only the electron in the 
final state is detectable), consideration of backgrounds 
will be crucial. 

Recent calculations of the neutrino-deuteron cross sec- 
tions are given in Refs. |3^|3^; in this energy range, the 
Ve cross section is ~ 65% larger than for Pg, and a factor 
~ 2 less than the inverse beta cross section on free pro- 
tons. The threshold for the i^e + d interaction is 1.4 MeV, 
less than the 4.0 MeV for Pe + d or 1.8 MeV for i>e + P- 
In the Ve channel, while the maximum electron (total — 
kinetic plus mass) energy is — E^ — 0.9 MeV, the peak 
in the differential cross section is lower [s^, , and so 
we use Ee — E^ ~ 2.5 MeV to assign energies. For the 
present purposes, a delta function with these kinemat- 
ics is an adequate approximation to the full differential 
cross section in the relevant energy range, and makes 
interpretation of the results more straightforward. The 
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total event rate is 



dE, 



dE^dE, 



(4) 



where N^i = 6 x 10 is the number of deuterons per 
kton of D2O, and T is the exposure time. In making our 
estimate of the SNO DSNB sensitivity, we assume that 
the total exposure will be 5 kton-yr at full efhciency. We 
emphasize that in order to more accurately determine 
the true sensitivity, a comprehensive study by the SNO 
collaboration is needed. 

In Fig. El we show our estimates of the DSNB Ve detec- 
tion spectra for three chosen temperatures. SNO is less 
sensitive than SK, and thus can only detect the DSNB if 
the flux is larger than expected. We have thus renormal- 
ized our standard predictions according to our detection 
criterion, justified in more detail below, of 1.6 expected 
signal events in 5 years with detected energy E^ between 
20 and 30 MeV. A detailed discussion of the backgrounds 
is needed, and we turn to that next. 
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B. Detector Backgrounds 

In Fig. 13 we show the renormalized DSNB spec- 
tra, along with our estimates of the most important 
background processes: solar neutrinos and sub-Cerenkov 
muons produced by atmospheric neutrinos. As we will 
show, the possible SNO sensitivity to DSNB Ve is sev- 
eral times worse than the existing SK limit on DSNB Ve\ 
therefore, we ignore DSNB interactions as a possible 
background. 

Since the proposed detection channel \s -\- d 
+ p + p, below about 20 MeV solar neutrinos provide 
an overwhelming background, even taking into account 
the fact that neutrino mixing reduces the flux of Ue- The 
solar neutrino fluxes are shown in Fig. ^ At the 
highest energies, the dominant background is from the 
solar hep reaction ('^He + p ^ ''He -I- e+ -I- i/^) 35]; even 
though its flux is ~ lO'' times smaller than the solar ^B 
beta decay ^3^] flux, its cndpoint is about 19 MeV, about 
3 MeV higher than for ^B. Since we only consider ener- 
gies above 20 MeV, this remains true even when energy 
resolution is taken into account, as we do. There is 
no possibility of significantly reducing this background, 
since it is the same detection reaction as the signal, and 
the angular distribution is only weakly backward (see 
Fig. 3 of Ref. [l^). Therefore we neglect consideration of 
all other backgrounds, and the precise details of the sig- 
nals, below 20 MeV. Additionally, we neglect the diffuse 
i>e flux made by all of the stars in the Universe before 
they end their lives (some as supernovae), since this is 
buried by the solar i>e flux 37]. 

In the SK DSNB search, the most significant 
background at high energies is from the electrons and 
positrons produced by sub-Cerenkov muons decaying at 
rest. These muons are invisible since they are produced 
inside the detector fiducial volume with non-relativistic 



FIG. 2: The DSNB signal spectra in SNO for a 5 kton- 
yr exposure, each renormalized to give an expectation of 
1.6 signal events in the electron energy range of 20-30 MeV 
(shaded). The renormalization factors for each temperature 
are given in Table |I] Below 20 MeV, solar neutrino interac- 
tions (the hep flux is shown with the dashed line) and other 
backgrounds are overwhelming. Above 30 MeV, the sub- 
Cerenkov (invisible) muon and other atmospheric neutrino 
backgrounds are too large. The energy resolution of SNO has 
been taken into account. In the irrelevant region below 20 
MeV, the signal curves are subject to some approximations 
in the calculation, and should not be taken too literally. 



initial energies, and they quickly stop; their decay spec- 
trum is simply the well-known Michel spectrum. Since 
SK has measured the rate of these decays, and their spec- 
trum is fixed, this greatly reduces the dependence on at- 
mospheric neutrino flux predictions. The invisible muons 
are produced by low-energy atmospheric i'^ and P^, pri- 
marily interacting with bound neutrons and protons in 
oxygen nuclei. After correcting for detector efficiency, 
the measured rate in SK is about 60 events per year in- 
side the 22.5 kton fiducial volume 0. Taking into ac- 
count vacuum neutrino oscillations of into i^r (and 
their antiparticles), this is in good agreement with ex- 
pectations iJi]. There are also background events due to 
very low-energy atmospheric and De, also primarily in- 
teracting with boimd nucleons. Since their rate is much 
smaller than the rate due to decaying invisible muons, 
especially below 40 MeV 0|, we do not consider them 
further. 

In the energy range 20-30 MeV, the efficiency- 
corrected rate of invisible muon decays in SK is about 
0.5/kton/yr. To estimate the rate of invisible muon de- 
cays in SNO, we have to correct for the fact that SK 
and SNO are at different latitudes. In the relevant neu- 
trino energy range, near 200 MeV, Ref. [l^ suggests a 
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TABLE L For each assumed temperature, we list the total 
DSNB flux and the flux just in the energy range 22.5 — 32.5 
MeV, both in units of cm~^ s~^. The last column gives the 
factor by which the flux would need to be increased to meet 
our detection criteria. 

T <f) (aU E,) <f) (22.5 - 32.5 MeV) (/)scns./0prcd. 



4 20 0.2 32 

5 16 0.4 15 

6 14 0.6 9 

7 12 0.8 7 

8 10 0.9 6 



correction of 1.8, while Ref. suggests 1.5. To be con- 
servative, we estimate that the rate will be a factor of 
2 larger at the location of SNO than in SK. We assume 
that events with a final-state neutron can be tagged, and 
that about half of the events produced on oxygen can 
be tagged by a nuclear de-excitation gamma from the 
residual ^^O or ^'^N ji^l; together these lead to a reduc- 
tion factor in the background rate of about 0.5, canceling 
the assumed latitude correction. In practice, the invisible 
muon decay rate will be measured and the assumed final- 
state probabilities checked. Our estimated residual invis- 
ible muon decay background rate in SNO of 0.5/kton/yr 
in the range 20-30 MeV is conservative, and will likely 
be reduced by several corrections we have neglected, each 
at the ~ 10 — 20% level and favorable. These include: a 
smaller latitude correction would be more accurate, the 
number of oxygen nuclei per kton of heavy water is less 
than in light water, capture on oxygen will reduce 
the number of background invisible muon decays, and 
the neutrino-deuteron cross sections are less than for free 
nucleons. Also, the observing period in SNO is mostly 
after that used in the SK analysis 4] , and closer to solar 
maximum, when atmospheric neutrino fluxes are smaller, 
particularly at high latitudes and low energies |4i| . 



C. SNO Sensitivity to the Electron Neutrino Flux 

To estimate the SNO sensitivity to DSNB Vf. in the 
energy interval 20-30 MeV, we consider two primary re- 
quirements. First, that in the assumed 5 kton-yr expo- 
sure, there should be at least one signal event. Second, 
that the number of signal events Ns should be larger 
than the square root of the number of background events 
Nb for the same exposure (after our assumed cuts). 
Since we estimate Nb = 2.5, the second requirement 
Ng > ^/Nb = 1.6 is hardly different from the first that 
Ng > 1. Our treatment is crude, but it is not yet possible 
to be more precise. The SNO backgrounds rates, which 
we have tried to estimate conservatively (see above) , have 
not yet been published. Further, since both Ng, Nb ^ 1, 
the statistical fluctuations in each will also be ~ 1, and 
what numbers of events happen to occur will affect the 



final conclusions. In addition, the energies at which they 
occur will matter - a single event at 21 MeV would be 
interpreted differently than a single event at 29 MeV. 
Since the statistical errors are so large, we can neglect 
all other uncertainties. In trying to motivate a complete 
and sophisticated analysis of the SNO data, we are sim- 
ply trying to establish the likely scale of the sensitivity 
and show why this would be an interesting and impor- 
tant result. The full analysis will have to be done by the 
SNO Collaboration. 

Since the sensitivity is above the level of standard pre- 
dictions, we consider how much larger the flux would have 
to be to yield a detection. This factor must be at least 

(?!>sens. max(l,V^) ,r\ 

0prcd. Ns 

where inside this equation, Ns is the number of signal 
events predicted in the standard case. Since Ns, Nb are 
so small, we neglect the possible spectral differences be- 
tween signal and background in the range 20-30 MeV, 
and only consider the counts, so that the required Ns 
does not depend on the assumed supernova neutrino tem- 
perature. Since the predicted flux in this range does de- 
pend on temperature, so does the required renormaliza- 
tion in order to have a detection. These results are given 
in Tabled as a function of the assumed temperature. For 
low temperatures, large fluxes would be required. On the 
other hand, for larger temperatures, the SNO sensitivity 
is reasonably close to standard predictions. We do not 
consider neutrino mixing for three reasons. First, the 
low statistics. Second, the fact that over this narrow en- 
ergy interval, a composite spectrum could be considered 
to be dominated by a largest temperature contributing 
to the spectrum. Third, since the estimated sensitivity 
is larger than standard predictions, any discovery would 
mean that our understanding of supernovae was missing 
an effect more significant than neutrino mixing. 

The estimated sensitivity can also be approximately 
characterized by the required flux in this interval, in- 
stead of the number of events. Multiplying the flux, the 
number of deuterons, the cross section at 25 MeV, and 
the assumed 5-year exposure, we find (6 x 10'^^) (27 x 
lO^'*^ cm^)(5 yrs) ~ 1.6 events, so that the required flux 
is 6 cm~^ s~^. Ignoring the slight temperature depen- 
dence, we show this as a constant in Fig. O 



IV. ELECTRON NEUTRINOS FROM SN 1987A 

The neutrino flux from SN 1987A was observed with 
two water Cercnkov detectors, Kamiokande-II (Kam-II) 
and 1MB 42, 43, 44, 45j. (And likely also with the much 
smaller Baksan scintillator detector, which had no direc- 
tional sensitivity ^3.) The most detectable component 
is the J>e flux, due to the largest cross section being the 
inverse beta reaction S'e+P — > + n, in which the outgo- 
ing positron angular distribution is nearly isotropic [231 ■ 
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Two features of the measured angular data remain 
mysterious. One is that the first event in Kam-II was 
directed forward. The other is that the majority of the 
12 Kam-II events, and all but one of the 8 1MB events, 
were in the forward hemisphere. Assuming that the yield 
is dominated by inverse beta events, both features, and 
especially the latter, are unlikely. These features suggest 
a larger than expected contribution from v + e^ — > i^ + e^ 
scattering, in which the electron closely follows the neu- 
trino direction. If so, then the most likely possibility is 
that this was caused by an enhanced flux of Uf,; the 
flux is constrained by the inverse beta yield, and the other 
flavors have ~ 6 times smaller cross sections on electrons. 

Taking the measured data at face value, a few inter- 
pretations of the noted features are possible. One, that 
the z/e flux in the early neutronization burst phase was 
much larger than expected. This could explain the first 
forward event, but not the time-averaged angular dis- 
tribution. Two, that the Ve flux in the longer thermal 
phase was much larger than expected. This could help 
explain the time-averaged angular distribution, and if the 
flux were large enough, also make it likely that one of 
the + — > i^e + e~ events would happen to be the 
first event detected. Third, that there were statistical 
fluctuations, either upward fluctuations of the small ex- 
pected yields of i^e + e~ ^ i^e + e~ scattering, or forward 
fluctuations of the near-isotropic angular distribution of 
De + P ^ e+ -t- n. With the SN 1987A data alone, it 
not possible to decide among these interpretations. In 
the context of standard supernova models, the first two 
interpretations are disfavored. What we consider here is 
that the required enhancements to the fluxes in the 
first two interpretations would similarly increase the ex- 
pectations for the DSNB I'e hux, and that this should be 
testable with SNO. 

At Kam-II the first event was perfectly forward, with 
an angle 18 ± 18 degrees relative to the neutrino direc- 
tion. During the neutronization phase, a prompt Vf, burst 
is created from electron captures on newly-liberated pro- 
tons, e- + i^e+n 47, ligi^lsj. The fact that 
the first event was both forward and first is suggestive 
of its being caused by neutrino-electron scattering [s^ . 
whatever the difficulties of explaining this with standard 
supernova models 0, E^- (The relative timing of 
the Kam-II and 1MB detectors is unknown, and we are 
making the common assumption that the first event in 
Kam-II was the first event overall.) The total energy re- 
leased during this burst is expected to be about an order 
of magnitude less than the Ve energy released during the 
thermal phase. For a total i^e burst energy of 4.5 x 10^^ 
erg j47|, the expected number of i^e scatterings at Kam- 
II is ~ 0.01 events, compared to ~ 0.1 for the thermal 
phase. These estimates neglect neutrino mixing, which 
can be especially relevant for the burst phase, since the 
transformation of V f, in to other flavors occurs without the 
reverse process [Sfl Isj, |53 . 

An enhanced rate of Ve + s~ ~> i^e + £~ during the ther- 
mal phase could help explain why the time-averaged an- 
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FIG. 3: The predicted DSNB flux as a function of the z/e 
temperature. The approximate Mont Blanc limit and pro- 
jected SNO sensitivity are shown with solid lines. The dotted 
lines approximately indicate the enhanced DSNB i/^ fluxes ex- 
pected if the apparent indications of z/e + e~ ^ z/e + e~ events 
in the SN 1987A data are interpreted as a probable outcome. 
The upper dotted line corresponds to assuming that the first 
Kam-II event was due to Ve from the neutronization burst 
phase. The lower dotted line corresponds to assuming in- 
stead that it arose from z/g from the thermal phase (and that 
in total ~ 3 of the events in Kam-II were due to z/g). 



gular distributions were so forward |56l l57| . A recent re- 
analysis determines that a few of the Kam-II events, and 
at most one of the 1MB events, are consistent with elastic 
scattering events jsH. With or without these events, the 
angular distributions of Kam-II and 1MB are not con- 
sistent with each other (or with expectations), when one 
considers both the mean and the variance of cos 9 '25', '5^. 
However, as a caution about the small-number statis- 
tics, note that had even a single additional event been 
detected in a backward direction in 1MB, then the dis- 
agreements would have been less severe [25| . 

Some potential scenarios for increased Ue emission have 
been considered in Refs. boost in the 

z^e yield relative to can be characterized by the factor 



(6) 



Here N is the total number of events observed at Kam-II 
or 1MB (12 or 8), n is the number of Ve + e~ ^ i^e + e~ 
scatterings, and (tr) is the thermally averaged cross- 
section, taking into account the thermal distribution and 
the efficiency of the detector ■ In deriving Eq. © , we 
have assumed that there is equal total energy emitted in 
each neutrino flavor. Since during the burst phase the 
total energy emitted in i/g is approximately an order of 
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magnitude less than by in the thermal phase, f^^ in- 
creases by i<;*^e''™aZ^^&Mrst ^jjgj^ considering a probable 

scenario of one event during the burst phase. 

For an illustration for a typical value of fi^^ for the 
thermal phase we consider n — i and Tp_, = 3 MeV for 
Kam-II. We choose n = 3 since this is the smallest n 
such that there is a reasonable chance of having an elec- 
tron scattering event be first, and also so that the time- 
integrated angular distribution is affected at a significant 
level. A typical value is /^^ ~ 15 over the range of T^^ . In 
Fig. O we show what this increase corresponds to in the 
DSNB flux, and how it compares to the SNO sensitivity. 
(If the z/e is assumed to come from the neutronization 
burst phase, the required boost is even larger, at least 
~ 100.) At high assumed temperatures, the suggested 
Ve flux enhancements in the thermal phase would lead 
to too many charged-current events on ^^O in the SN 
1987A data, which would be an independent way to con- 
strain these scenarios. For example, for standard fluxes, 
this yield would become a few events when T = 8 MeV; 
in addition, the angular distribution favors backward an- 
gles ^|. Also, for large temperatures and fluxes, the 
number of events expected in the Homestake detector 
might become too large |6^ . 

There also remains the possibility that the first Kam-II 
event resulted from Ve+P +n, and just happened to 
be emitted in the forward direction. For example, there 
is a ^ 5% chance that any 9^ event was emitted within a 
given cone of half-angle 25 degrees, comparable to the un- 
certainty in the Kam-II angular resolution. To have the 
first event be forward thus has a 5% probability, and 
to explain the time-averaged angular distribution would 
require additional statistical fluctuations. By searching 
for the DSNB flux, SNO will be the first experiment 
able to revisit the electron-neutrino component of the SN 
1987A data and determine if the observed angular distri- 
bution was indeed a probable outcome. 



V. CONCLUSIONS 

We have shown that SNO should have sensitivity to 
a Diffuse Supernova Neutrino Background (DSNB) I'e 
flux as low as ~ 6 cm~^ s~^, which is about three or- 
ders of magnitude smaller than the current limit from 
the Mont Blanc experiment In order to best avoid 
detector backgrounds, this is the flux just in the energy 
range 22.5-32.5 MeV; our corresponding theoretical pre- 
dictions, given in Table HI are ^ 1 cm~^ s~^. Our results 
depend on our estimates of the relevant solar and atmo- 
spheric neutrino backgrounds, and how the latter might 
be reduced. To the extent possible, we based our es- 
timates on measured rates, especially the sub-Cerenkov 
muon decay rate from SK. Given the very small projected 
statistics, a more precise estimate of the sensitivity is not 
yet possible. A full analysis by the SNO Collaboration, 
using the measured data, is strongly encouraged. 

Our estimate for the SNO sensitivity to DSNB i/g may 



seem surprisingly good, given that SNO is much smaller 
than SK. However, it can be conflrmcd by adjusting the 
existing SK DSNB limit of 1.2 cm^^ s^^ for the dif- 
ference in exposure, ignoring the small differences in the 
considered energy ranges. The SK exposure was 22.5 
kton for 4.1 years at about 50% efficiency in the rele- 
vant energy range, so about 45 kton-yr at full efficiency, 
or about 9 times more than the assumed SNO expo- 
sure. Since the SK limit arises from a background-limited 
search, the scaling to 5 kton-yr in a light water detector 
would worsen the limit by about -^79 = 3. Taking into 
account the smaller signal cross section in a heavy water 
detector (and assuming the same background rate) gives 
an additional penalty of a factor ~ 2, so that in the end 
the SNO sensitivity should be about a factor of 6 worse 
than in SK, close to what we deduced. (Had SK been 
rate-limited, the estimated SNO sensitivity would have 
instead been about 18 times worse.) Note that if the 
true SNO exposure is less than the 5 years of full effi- 
ciency that we assumed, the flux sensitivity only changes 
by the square root of the ratio of exposures. 

An interesting related point is that once the Kam- 
LAND and SNO detectors report on multiple kton-yr 
exposures, they should be able to set DSNB limits 
only a factor of a few to several weaker than what SK 
has published so far. This would be significant princi- 
pally because KamLAND and SNO should have lower 
backgrounds and thresholds; the combination of these 
limits with the stronger limits at higher energies in SK 
might lead to irnproved limits on the supernova rates at 
higher redshiftJ^llS. The KamLAND and SNO hm- 
its so far |3lL Im | are significantly worse than we are 
estimating here because they were based on the small 
exposures of early data and considered narrower energy 
ranges. In addition, the SNO analysis was based on the 
pure D2O phase, during which the neutron detection ef- 
ficiency was much lower than in the next two phases. 
The sensitivity improves linearly with the exposure be- 
fore the background-dominated regime is reached, after 
which it improves only with the square root. In terms 
of absolute sensitivity, the signal requirements demand 
an exposure of at least several kton-yr, beyond which we 
estimate that SNO will be background-dominated. (On 
a related issue, note that while the present SK DSNB 
sensitivity is background-limited, if SK is enhanced by 
the addition of gadolinium, then the absolute sensitiv- 
ity will be greatly improved, and is projected to become 
rate-limited 9].) 

Studying the DSNB with any Ve channel other than 
i^e + d requires a large cross section, large detector mass, 
and favorable kinematics. In the future, this may be 
possible with large liquid argon detectors, as long empha- 
sized by Cline 65] . For example, the estimated DSNB 
i^e sensitivity after a 15 ktoii-yr exposure of the ICARUS 
detector is ~ 1.6 cm~^ |66l |. which would be very in- 
teresting as a complement to the SK DSNB sensitiv- 
ity. This estimate assumes that the only backgrounds in 
ICARUS are those due to solar and atmospheric i^e ■ Very 
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likely, other backgrounds will have to be considered as 
well. For example, it is assumed that the low-momentum 
(sub-Cerenkov) muons, which are a serious background 
in water-based detectors, could be completely rejected in 
argon-based detectors |6a,|6y|- In water-based detectors, 
the muons with momenta below ~ 100 MeV are invisible. 
In a liquid argon detector with an assumed kinetic energy 
threshold of 5 MeV, muons with momenta less than ~ 30 
MeV would also be invisible, and their decays would pro- 
duce a relevant background. Also needed are studies of 
backgrounds induced by fast neutrons entering the de- 
tector, as well as those arising from the quenched light 
output of energetic charged particles. Nevertheless, this 
technique appears promising, and may have the potential 
to go significantly beyond the sensitivity of SNO. 

Besides waiting for another Milky Way supernova, in 
principle there is another method to probe the t^e emis- 
sion. The flux from past Milky Way supernovae would 
have transformed terrestrial molybdenum isotopes into 
unstable technetium isotopes (half-lives 2.6 and 4.2 mil- 
lion years) , and the elemental ratios in deep molybdenum 
deposits should reveal the received v,, fluence over a com- 
parable time period [b^]. To the best of our knowledge, 
the realistic sensitivity of this technique is unknown, and 
no experiment is planned. 

In addition to improving the present DSNB v^, sensi- 
tivity by about three orders of magnitude, a full analysis 
by SNO can test whether the large flux suggested by 



the SN 1987A data was probable, or was simply a sta- 
tistical fluctuation. If the first event in Kamiokande-II 
was really Ve + e~ ^ Ve + e~ scattering, and was a prob- 
able outcome, then the DSNB Vf. flux should be greatly 
enhanced. If this first event is assumed to be from the 
neutronization phase, then the required enhancement is 
a factor at least ~ 100. If this first event is assumed to 
be from the thermal phase, the required enhancement is 
a factor at least ~ 15; this would also mean that there 
should be at least a few other events in the data, which 
would make the angular distributions more forward, also 
as suggested by the data. Both scenarios are probably 
unrealistic, and our estimates are crude. However, our 
point is that at present, only SNO has the sensitivity to 
directly constrain them, and for the first time be able to 
look back and shed light on the mysteries still lingering 
from SN 1987A. 
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